SYSTEMATIC SEARCH FOR SHORT-TRANSIENTS AND PULSATION EVENTS FROM 

INTEGRAL SURVEY DATA 



o 
o 

(N 



Ken Ebisawa, Peter Kretschmar, Nami Mowlavi, Ada Paizis, Nicolas Produit, Simon Shaw^, 
Sandro Mereghetti, Diego Gotz'^, Stefan Larsson'^, Niels Joergen Westergaard^, Sami Maisala^, 

and Riidiger Staubert^ 

^INTEGRAL Science Data Centre, Chemin d'Ecogia 16, 1290 Versoix, Switzerland 

^CNR-IASF, via Bassini 15, 1-20133 Milano, Italy 

■^Stockholm Observatory, AlbaNova, SE-106 91 Stockholm, Sweden 

^Danish Space Research Institute, Juliane Maries Vej 30, DK-2100 Copenhagen 0, Denmark 

^Observatory, P.O. Box 14, FIN-00014 University of Helsinki, Finland 
^Institut fiir Astronomic und Astrophysik, Astronomic, Sand 1, D-72076 Tiibingen, Germany 



(N 



> 

\o 

o 
^■ 

o; 
^. 

:^: 

o 



X 



ABSTRACT 



The imaging instruments on board INTEGRAL have 
wide fields of view and high time resolution. There- 
fore, they are ideal instruments to search for pulsat- 
ing sources and/or transient events. We are system- 
atically searching for pulsations and transient events 
from known and serendipitous sources in the Galactic 
Plane Scan (GPS) and Galactic Center Deep Expo- 
sure (GCDE) core program data. We analyze the 
standard pipe-line data using ISDC Off-line Science 
Analysis (OSA) system, so that our results are re- 
producible by general guest users. In this paper, we 
describe our system and report preliminary results 
for the first year of operation. 
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1. SCIENTIFIC OBJECTIVE 



transient events, which are likely to be smeared if 
integrated for the entire ScW. Also, QLA does not 
search for pulsations, and fast pulsars down to ~msec 
periodicity may not be detectable unless event arrival 
time analysis is performed; such an analysis is not 
considered in the standard analysis pipe-line. 

The INTEGRAL Burst Alert System (IB AS; 
Mereghetti 2003) produces real-time results of bright 
transient events; not only Gamma-ray bursts, but 
also other transient events at various timescales from 
milliseconds to minutes. However, IBAS directly ac- 
cesses raw telemetry, and is not able to read the stan- 
dard processed data. Also IBAS is not available to 
public. 

To supplement the QLA, standard analysis and 
IBAS, we are developing a system to systematically 
search for transient events and pulsations from core 
program GPS and GCDE data. Our system will be 
publicly available as a part of ISDC OSA. In this pa- 
per, we describe the technical aspects of our project, 
and demonstrate the feasibility. 



The wide fields of view and the dithering strategy 
of INTEGRAL observations enable us to carry out a 
systematic survey of a large area of the sky. The IN- 
TEGRAL Science Data Center (ISDC; Courvoisier 
ct al. 2003), monitors the observations of the instru- 
ments with the Quick Look Analysis (QLA) software. 
This allows us to systematically search for bright 
transient events, as well as monitoring the primary 
targets of the observation, within a few hours of re- 
ceiving the telemetry (Shaw et al. 2004). 

However, the QLA carries out analysis only on the 
Science Window (ScW) basis. Each ScW lasts about 
30 minutes, and analysis at shorter timescales, which 
would take too much CPU time and disk space, is 
not feasible with QLA. Therefore, we may miss short 



2. ANALYSIS METHODS 



2.1. Transient events from known sources 



We routinely run the standard pipeline analysis on 
the GPS and GCDE data. For the known and de- 
tected sources, we make ISGRI and JEM-X light 
curves with 100 second time-bin^. 

After the JEM-X and ISGRI light curves are made, 
we carry out statistical tests to search for variable 



^We use ISDC cxccutables "iijight" and "j_srcjc" to create 
ISGRI and JEMX source light curves, respectively. 
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Figure 1. Example of JEM-X light curves of variable sources, which show significant variations within a single 
Science Window. Bin width is 100 seconds, and the light curves are shown for each Science Window. 
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Figure 2. Left: ISGRI Vela X-1 light curve from TJD = 12970 to 13006 (where Truncated Julian Date = MJD 
- 40000). Bin-width is 512 sec, and the energy range is from 20 keV to 60 keV. Right: Expansion of the large 
flare at around TJD =12971. Bin-width = 128 sec. 



events. In Figure 1, we show examples of JEM-X 
light curves with significant variations within one 
ScW. In this case, the search for variation is made 
for each ScW. Namely, we find only variations shorter 
than the length of a ScW. 

Since each Sew is a single, stable pointing, these 
short-term variations are free from systematic effects 
that may be introduced by the observation dithering 
pattern. 

In Figure 2, we show an example of the ISGRI Vela 
X-1 light curve. A prominent flare is seen at around 
TJD (Truncated Juhan Date) = 12971 (Krivonos et 
al. 1994; Staubert et al. 2004). 



2.2. New transient search 



While IBAS systematically searches for strong 
gamma-ray bursts with IBIS, soft X-ray transients 
below ^15 keV may be expected only by JEMX. 
Also, IBAS may not trigger weak gamma-ray bursts, 
which may be found in the off-line analysis. 

To search for new short-term transients, which are 
not in the reference catalog, we need to carry out 
image deconvolution and do a source search for a 
much shorter duration than the ScW. Current OSA 
(v. 3.1) is not really suitable for such a purpose, since 
a user Good Time Interval (GTI) has to be specified 
for each Observation Group before data correction 
and instrumental GTI calculation are made. 

We adopt a new, more efficient data processing ffow, 
which will be standard in the INTEGRAL archives 
and OSA 4.0 (Figure 3). In the new scheme, after 
the standard RAW, PRP, COR and GTI process, 
a merged event list is created for each ScW. Each 
event has physical arrival time (TIME column), as 
well as START and STOP times of GTI. The GTI 
extensions are included in the same file, therefore 
generic tools (such as xronos ftools, xselect) can read 
these event files. 

After the event list is made for a ScW (which is typ- 
ically 30 minutes or so), we define many short user 
GTIs, each 1 minute or so. By applying the short 
user GTI, we run OSA and create sky images for 
each ~ 1 minute, and carry out a source search. Note 
that in this method the user GTI is specified after the 
Correction and instrumental GTI calculation step, so 
that we do not have to duplicate these processes for 
each user GTI. 

We have already adapted OSA so that it runs on 
the merged event file, and have started to systemati- 
cally search for new transient events in the GPS and 
GCDE data. We are evaluating sensitivity of ISGRI 
and JEMX for short transient events. Preliminary 
results suggest we may detect ~ 50 mCrab transient 
source for one minute interval with both JEMX (5 - 
30 keV) and ISGRI (20 - 60 keV). 



2.3. Short pulsation search 

To search for fast (e.g. ^ millisecond) pulsations, we 
need to make event arrival time analysis without im- 
age deconvolution, since it is not practical to carry 
out deconvolution for such a short period. The ISDC 
tool "barycent" (since version 2.1) has been modified 
so that it can perform the barycentric correction on 
the event lists. To optimize the search for pulsations, 
we select events based on the Pixel Illumination Fac- 
tor (PIF), which is the fractional area of each ISGRI 
pixel exposed to a source. The PIF takes values from 
to 1, such that PIF=0 means this pixel is not sen- 
sitive to that source, and PIF=1 means this pixel is 
fully exposed to the source. 

A standard and stand-alone tool ii_pif_build is 
used to calculate the ISGRI PIF for a given source. 
Once the PIF map is made, it is straightforward to 
filter event lists based on the PIF threshold (Figure 
3). In Figure 4, we show an example of the PIF and 
PIF-based event filtering. By applying the PIF filter, 
we can reduce the background and optimize the sig- 
nal to noise ratio. From experience, we choose PIF 
> 0.5 as the event selection criterion. A similar PIF 
creation tool is expected for JEMX too. 

In Figure 5, we demonstrate the effect of PIF selec- 
tion. In this example, we analyze ISGRI Crab pul- 
sar data from ScW 010200380060. The events are 
folded with the apparent pulse-period (barycentric 
correction is not made) at 33.550 msec. We can see 
that, after the PIF selection, the background is sub- 
stantially reduced and the pulse-shape is much more 
prominent. Note that the fmaskf ilt ftool is used to 
carry out the ISGRI PIF selection. 

We have confirmed our event timing analysis method 
by detecting pulsations from known pulsars such as 
Cen X-3. By combining JEMX and ISGRI, we will 
be able to study known pulsars and search for new 
pulsars in the wide energy range from ~ 5 keV to 
- 200 keV. 
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Figure 3. Analysis flow for our short transient search and pulsar search program. 
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Figure 4- Lefl: ISGRI pixel Illumination Factor (PIF). Middle: Original detector image. Right: Detector image 
with pixels whose PIF values exceed 0.5. The data is from Science Window 010200380060 (Crab observation). 
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Figure 5. Folded pulse-profile of the Crab pulsar from ISCRI ScW 010200380060. Black is made from the 
original event file using all the events, and red is from the events whose PIF values exceed 0. 5. Each profile is 
normalized so that the pulse- fractions relative to the average counting rates are plotted. 



